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ABSTRACT 

We present redshift evolution of galaxy effective radius obtained from the Hubble Space Telescope 
{HST) samples of ~ 190, 000 galaxies at z = 0 — 10. Our HST samples consist of 176,152 photo-z 
galaxies at z = 0 — 6 from the 3D-HST+CANDELS catalogue and 10,454 Lyman break galaxies 
(LBGs) at z = 4 — 10 identified in CANDELS, HUDF09/12, and HFF parallel fields, providing the 
largest data set to date for galaxy size evolution studies. We derive with the same technique over 
the wide-redshift range of z = 0 — 10, evaluating the optical-to-UV morphological A-correction and 
the selection bias of photo-z galaxies-|-LBGs as well as the cosmological surface brightness dimming 
effect. We find that Vg values at a given luminosity significantly decrease towards high-z, regardless 
of statistics choices (e.g. rg oc (1 -I- for median). For star-forming galaxies, there is no 

evolution of the power-law slope of the size-luminosity relation and the median Sersic index (n ^ 1.5). 
Moreover, the rg-distribution is well represented by log-normal functions whose standard deviation 
cinre does not show significant evolution within the range of ainr^ ~ 0.45 — 0.75. We calculate the 
stellar-to-halo size ratio from our measurements and the dark-matter halo masses estimated from 
the abundance matching study, and obtain a nearly constant value of rg/rvir = 1.0 —3.5% at z = 0 — 8. 
The combination of the rg-distribution shape-l-standard deviation, the constant Tg/rvir, and n ^ 1.5 
suggests a picture that typical high-z star-forming galaxies have disk-like stellar components in a sense 
of dynamics and morphology over cosmic time of z ^ 0 — 6. If high-z star-forming galaxies are truly 
dominated by disks, the rg/rvir value and the disk formation model indicate that the specific angular 
momentum of the disk normalized by the host halo is jd/nid — 0.5 — 1. These are statistical results 
for galaxies’ major stellar components, and the detailed study of clumpy sub-components is presented 
in the paper H. 

Subject headings: cosmology: observations — early universe — galaxies: formation — galaxies: high- 
redshift 


1. INTRODUGTION 

Galaxy sizes offer a variety of invaluable insights into 
the galaxy formation and evolution. A slope of size- 
stellar mass (or luminosity) relation, a size growth rate, 
and a size distribution are key quantities for understand¬ 
ing developments of galaxy morphology and properties of 
host dark matter (DM) halos. 

Studies of high-z galaxy sizes show substantial pro¬ 
gresses by observations of the Hubble Space Telescope 
(HST) which is capable of high spatial resolution imag¬ 
ing. Galaxy sizes defined by the effective radius, rg, have 
extensively been measured with Advanced Gamera for 
Survey (ACS) and Wide Field Camera 3/IR-channel on¬ 
board HiST for massive galaxies at 0 < z < 3 (e.g., van 
der Wei et al. 2014) and z > 3 — 4 Lyman break galaxies 
(LBGs) selected in the dropout technique (Steidel et al. 
1999) (e.g., Trujillo et al. 2006; Toft et al. 2007, 2009; 
Dahlen et al. 2007; Grazian et al. 2012; Huang et al. 2013; 
McLure et al. 2013). However, these studies, particularly 
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at high-z, do not reach agreement on the size growth 
rate. Oesch et al. (2010) have reported that the average 
size evolves according to roughly rg oc (1 -|- z)“^ based 
on a z ^ 7 LEG sample in an early-epoch data of their 
HST survey (see also e.g., Bouwens et al. 2004; Holwerda 
et al. 2014). On the other hand, Hathi et al. (2008b) have 
argued that the average size scales asrg oc (l-|-z)“^-® 
using LBGs at z ~ 2 — 6 (see also e.g., Ferguson et al. 
2004). Some studies have provided results of a growth 
rate falling between these two growth rates (e.g., Mosleh 
et al. 2012, 2013; Ono et al. 2013). Moreover, Gurtis- 
Lake et al. (2014) have suggested no significant evolution 
of typical galaxy sizes, if one does not use average but 
modal values of size distribution for representative radii 
at a given redshift. These discrepancies in the evolution¬ 
ary trend would be attributed to small HST samples at 
z > 3 —4 and/or potential biases caused by heterogenous 
samples and measurements taken from the literature. 

The two size growth rates of Cg oc (1 -|- z)“^ ® and 
Tg oc (1 -I- z)~^ correspond to the cases of a fixed virial 
mass and circular velocity of DM halos, respectively, if 
the stellar-to-halo size ratio (SHSR) is constant over the 
redshift range. Assuming the constant SHSR, a number 
of studies discuss the evolution of host DM halos with 
the size growth rates (e.g., Ferguson et al. 2004; Hathi 
et al. 2008a). However, the evolution of SHSR is not 
well understood. Recently, SHSRs have been estimated 
observationally with results of abundance matching tech- 
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TABLE 1 

Limiting Magnitudes of HST Images for Size Analyses 


Field 

(1) 

Veoe 

(2) 

7814/2850 

(3) 

15cr (5(t) Depth 
Ji25 Hieo 

(4) (5) 

Pb98TL05 4l25^^160^ 
(6) 

Ji25Hieo^ 

(7) 

HUDF09-I-12 

29.3 (30.5) 

28.1 (29.3) 

28.7 (29.9) 

28.7 (29.9) 

29.4 (30.6) 

29.2 (30.4) 

HUDF09-P1 

29.4 (30.6) 

28.7 (29.9) 

28.0 (29.2) 

27.8 (29.0) 

28.4 (29.6) 

28.2 (29.4) 

HUDF09-P2 

28.2 (29.4) 

27.5 (28.7) 

28.2 (29.4) 

28.0 (29.2) 

28.6 (29.8) 

28.5 (29.7) 

GOODS-S Deep 

27.6 (28.8) 

27.6 (28.8) 

27.2 (28.4) 

27.1 (28.3) 

27.8 (29.0) 

27.5 (28.7) 

GOODS-S Wide 

27.6 (28.8) 

27.2 (28.4) 

26.6 (27.8) 

26.4 (27.6) 

27.1 (28.3) 

26.9 (28.1) 

GOODS-N Deep 

27.6 (28.8) 

29.3 (30.5) 

27.1 (28.3) 

26.9 (28.1) 

27.3 (28.5) 

27.4 (28.6) 

GOODS-N Wide 

27.5 (28.7) 

28.7 (29.9) 

26.4 (27.6) 

26.3 (27.5) 

26.9 (28.1) 

26.7 (27.9) 

UDS 

27.0 (28.2) 

27.0 (28.2) 

26.3 (27.5) 

26.4 (27.6) 

26.7 (27.9) 


COSMOS 

27.1 (28.3) 

26.8 (28.0) 

26.4 (27.6) 

26.4 (27.6) 

26.7 (27.9) 


AEGIS 

27.1 (28.3) 

26.6 (27.8) 

26.4 (27.6) 

26.5 (27.7) 

26.8 (28.0) 


HFF-Abell2744P 

27.9 (29.1) 

27.6 (28.8) 

27.6 (28.8) 

27.6 (28.8) 

28.1 (29.3) 

27.9 (29.1) 

HFF-MACS0416P 

27.7 (28.9) 

27.6 (28.8) 

27.9 (29.1) 

27.8 (29.0) 

28.3 (29.5) 

28.1 (29.3) 

PSFFWmr 

0708 

0709 

0(T2 

(FTS 

(FTB 

(FTB 


Note. — Columns: (1) Field. (2)-(7) Limiting magnitudes defined by a 15(T (Sct in parentheses) sky noise in a 
0^^35-dianieter aperture. 

^ Stacked image of YqqsFios Ji 25 F^i 60 -t>ands for LBGs at 2 ~ 4,5,6, and 7. The JH 140 image is also included for 
LBGs at 2 ~ 7 — 8 in the HUDF09+12 field. 

^ Stacked image of J 125 - and ffisg-bands for LBGs at 2 ~ 8. 

Typical size of PSF FWHMs. 


niques (e.g., Behroozi et al. 2010, 2013) for galaxies at 
z ~ 0 (Kravtsov 2013) and at z ~ 2—10 (Kawamata et al. 
2014). Kawamata et al. (2014) conclude that virtually 
constant value of SHSR, 3.3±0.i%, over the wide redshift 
range. Galaxy disk formation models of e.g., Fall (1983, 
2002); Barnes & Efstathiou (1987); Mo et al. (1998) pre¬ 
dict that galaxy disks acquire an angular momentum 
from its host DM halo trough tidal torques during the 
formation of these systems, leading to the proportional¬ 
ity between the two sizes. The SHSR values provide us 
with information about the DM spin parameter and the 
fraction of specific angular momentum transferred from 
DM halos to the central galaxy disks (e.g.. Mo et al. 
1998). 

Additionally, the size-stellar mass relation and the 
scatter of size distribution present independent evidence 
for the picture of galaxy disk formation (e.g.. Fall 1983, 
2002; Shen et al. 2003; Bullock et al. 2001). van der Wei 
et al. (2014) have revealed that the slope of size-stellar 
mass relation and the scatter do not significantly evolve 
0 < z < 3 in a systematic structural analysis for large 
samples of star-forming galaxies (SFGs) and quiescent 
galaxies (QGs) with a photometric redshift (photo-z). 
The constant values of these quantities strongly suggest 
that the sizes of SFGs are determined by their host DM 
halos. However, the controversial results of the slope and 
scatter evolution are obtained at z > 3 — 4 (e.g., Huang 
et al. 2013; Gurtis-Lake et al. 2014), probably due to 
large statistical uncertainties given by the small galaxy 
samples. An analysis with a large LBG sample would 
reveal the galaxy structure evolution up to z ^ 10 with 
no significant statistical uncertainties, and allow us to 
understand disk formation mechanisms, internal star for¬ 
mation, and morphological evolution over cosmic time. 

In this paper, we systematically investigate redshift 
evolution of galaxy sizes with an unprecedentedly large 
sample of 186,603 galaxies at z = 0 — 10 made from 
the HST deep data of extra-galactic legacy surveys. We 
assess effects of morphological A-correction, statistics 
choice, and sample selection bias with the galaxies at 
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Fig. 1.— Relation between UV magnitude and stellar mass for 
the SFGs at 2 0 — 6. The circles with error bars indicate SFGs 

at 2 = 0 — 1 (red), 2 = 1 — 2 (magenta), 2 — 3 (orange), 3 — 4 
(green), 4 — 5 (cyan), and 5 — 6 (blue). The colored solid lines 
present the best-fit relation at each 2 bin with the color 

coding same as the one of the circles.The black solid line denotes 
the M*-Muv relation for LBGs at 2 4 in Gonzalez et al. (2011). 

The top x-axis provides the corresponding UV luminosity in units 
of L*_^. The error bars denote the 16th and 84th percentiles of 
distribution. The vertical dashed lines denote thresholds of L\jy 
bins, 1, 0.3, 0.12, and 0.048 Luv/^ 2 _ 3 , f^^m left to right. 

z < 4, and then extend our systematic morphological 
measurements to 2 > 4. This paper has the following 
structure. In Section 2, we describe the details of our 
HST galaxy samples. Section 3 presents methods for es- 






















Evolution of Galaxy Size at z = 0 — 10 
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TABLE 2 

Number of Photo-2 Galaxies for our Size Measurements 


Field 

(1) 

2 = 0-1 
(2) 

2 = 1 - 2^= 

(3) 

AgALFIT / AsfG , QGs 

2 = 2-3'" 2 = 3- 4 

(4) (5) 

z = 4 — 5 
(6) 

z = 5 - 6 
(7) 

SFGs, r?P‘(4500 - SOOOA) 

HUDF09+12 

294 (397) 

368 (611) 

8 (157) 




HUDF09-P'‘ 

2168 (3467) 

2024 (4707) 

66 (1451) 




GOODS-S Deep 

1753 (2793) 

1402 (3847) 

37 (l296) 




GOODS-S Wide 

2790 (4724) 

2267 (6313) 

66 (25I8) 




GOODS-N Deep 

3270 (5106) 

1778 (5I68) 

71 (2094) 




GOODS-N Wide 

3903 (5939) 

2731 (6611) 

139 (2477) 




UDS 

5157 (9175) 

5433 (15771) 

209 (6094) 




AEGIS 

6441 (11074) 

5833 (13943) 

278 (6418) 




GOSMOS 

6856 (ll385) 

3594 (9754) 

179 (3915) 




-^totai(^) 

32632 (54060) 

25430 (66725) 

1053 (26420) 




-^total 



59115 (147205) 



SFGs, ryV(i5oo-3000A) 

HUDF09+12 


145 (611) 

79 (157) 

34 (69) 

19 (33) 

12 (26) 

HUDF09-P=" 


777 (4707) 

624 (1451) 

432 (936) 

160 (453) 

102 (177) 

GOODS-S Deep 


776 (3847) 

633 (l296) 

348 (696) 

101 (347) 

40 (702) 

GOODS-S Wide 


1297 (6313) 

1154 (2518) 

535 (1147) 

138 (487) 

44 (213) 

GOODS-N Deep 


1235 (5I68) 

784 (2094) 

389 (987) 

154 (446) 

66 (174) 

GOODS-N Wide 


1711 (6611) 

1114 (2477) 

412 (962) 

165 (5I6) 

47 (l67) 

UDS 


2730 (15771) 

1747 (6094) 

678 (2266) 

180 (7I6) 

52 (l76) 

AEGIS 


3158 (13943) 

2182 (6418) 

952 (2768) 

228 (873) 

84 (281) 

GOSMOS 


2413 (9754) 

1642 (3915) 

939 (2048) 

192 (765) 

54 (296) 

-^total(^) 


14242 (66725) 

9959 (26420) 

4719 (ll879) 

1337 (4636) 

501 (l796) 

-^total 



30765 (165517) 





QGs, 

(4500 - 8000A)‘" 




HUDF09+12 

323 (743) 

133 (458) 

2IM) 




HUDF09-P=‘ 

267 (637) 

113 (365) 

1 (72) 




GOODS-S Deep 

193 (447) 

85 (261) 

1 (46) 




GOODS-S Wide 

369 (895) 

115 (444) 

4 (81) 




GOODS-N Deep 

259 (623) 

86 (346) 

2 (98) 




GOODS-N Wide 

272 (744) 

110 (382) 

2 (147) 




UDS 

320 (1375) 

221 (933) 

7 (207) 




AEGIS 

387 (l369) 

270 (824) 

6 (272) 




GOSMOS 

890 (l738) 

170 (602) 

7 (127) 





3013 (7934) 

1190 (4250) 

31 (l076) 




-^total 



4234 (13260) 




Note. — Columns: (1) Field. (2)-(7) Number of the photo -2 galaxies that have S/N > 15 and reliable GALFIT 
outputs in each redshift range. The value in parentheses is the number of the photo -2 galaxies in the parent sampie. 


^ Total number of objects in the HUDF09-P1 and HUDF09-P2 fieids. 

^ The actual redshift range is 2 < 2 < 2.1 (1.2 < 2 < 2) for the measurement. 

^ The numbers of QGs with are not shown here due to the rarity at 2 > 2 — 3 and the UV faintness. 


timating galaxy sizes. In Section 4, we evaluate the mor¬ 
phological K corrections, statistics-choice dependences, 
and selection biases. We show the redshift evolution of 
size-relevant physical quantities in Section 5. Section 6 
discusses the implications for galaxy formation and evo¬ 
lution with results of our structural analyses. We sum¬ 
marize our findings in Section 7. 

Throughout this paper, we adopt the concordance 
cosmology with (fl^, ^A, h) = (0.3,0.7, 0.7), (Ko¬ 

matsu et al. 2011). All magnitudes are given in 
the AB system (Oke & Gunn 1983). We re¬ 
fer to the HST F606W, F775W, F814W, F850LP, 
F098M, F105W, F125W, F140W, and F160W filters as 

heoe! *775:-^814,-^850: ho98, Fio 5: <^125: di7i4o, ^nd i7i60: re¬ 
spectively. 

2. DATA AND SAMPLES 

We make use of the following two galaxy samples con¬ 
structed from the deep optical and near-infrared imag¬ 
ing data taken by ffST deep extra-galactic legacy sur¬ 


veys whose limiting magnitudes and PSF FWHM sizes 
are summarized in Table 1. In the last subsection, we 
explain the stellar masses of the sample galaxies. 


2.1. Sample of Photo-z Galaxies at z = 0 — 6 
in 3D-HST+CANDELS 

The hrst sample is made of 176,152 HST/WYCi-IR, 
detected galaxies with photometric redshifts (hereafter 
photo-z galaxies) at z = 0 — 6 taken from Skelton et al. 
(2014). These galaxies are identified in five Gosmic As¬ 
sembly Near-infrared Deep Extragalactic Legacy Survey 
(CANDELS) fields (Grogin et al. 2011; Koekemoer et al. 
2011), and detected in stacked images of Ji 25 ,JHi 4 o, 
and Hiqq bands of WFG3/IR, which yields roughly a 
stellar mass-limited sample. The photometric proper¬ 
ties and the results of spectral energy distribution (SED) 
htting for all the sources are summarized in Skelton 
et al. (2014). The HST images and catalogues are pub- 
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TABLE 3 

Number of LBGs for our Size Measurements 


Field 

(1) 

2 4 
(2) 

z b 
(3) 

Ngalfit /Nlbg 

2 ~ 6 2 ~ 7 

(4) (5) 

z 8 
(6) 

z 10 

(7) 

HUDF09+12 

160 (348) 

43 (130) 

26 (86) 

13 (50) 

9 (24) 

0 (2) 

HUDF09-P1 


41 (95) 

12 (30) 

2 (9) 

2 (7) 

0 (0) 

HUDF09-P2 


30 (90) 

8 (37) 

4 (23) 

0 (16) 

0 (0) 

GOODS-S Deep 

1046 (1872) 

292 (696) 

122 (311) 

55 (203) 

11 (57) 

1 (1) 

GOODS-S Wide 

294 (510) 

73 (142) 

20 (51) 

9 (31) 

3 (21) 

0 (0) 

GOODS-N Deep 

868 (1655) 

279 (630) 

48 (135) 

35 (ill) 

12 (28) 

1 (2) 

GOODS-N Wide 

522 (800) 

106 (222) 

25 (68) 

12 (231) 

3 (28) 

1 (1) 

UDS 


152 (310) 

39 (65) 

12 (25) 



AEGIS 


189 (381) 

47 (101) 

11 (28) 



GOSMOS 


209 (348) 

40 (80) 

11 (27) 



HFF-Abell2744P 


15 (37) 

12 (26) 

4 (7) 

2 (7) 


HFF-MACS0416P 


30 (134) 

23 (106) 

5 (18) 

4 (10) 


-^total {^) 

^total 

2890 (5185) 

1459 (3215) 

422 (1096) 

4993 (10454) 

173 (763) 

46 (l95) 

3 (6) 


Note. — Columns: (1) Field. (2)-(7) Number of the LBGs that have S/N > 15 and reliable 
GALFIT outputs in each redshift range. The value in parentheses is the number of LBGs in the parent 
sample. 


licly released at the 3D-HST website®. The catalogues 
include the spectroscopic redshifts on the basis of the 
HST/WFC3 G141 grism observation (Brammer et al. 
2012). We use galaxies whose physical quantities and 
photometric redshifts are well derived from SED fitting 
(specifically, sources with use_phot= 1 in the public cat¬ 
alogues). Tables 2 summarizes the number of galaxies at 
each redshift in the photo-z galaxy sample that we use. 
In this paper, we assume Salpeter (1955) initial mass 
function (IMF). To obtain the Salpeter IMF values of 
stellar masses (M*) and star formation rates (SFRs), we 
multiply the Chabrier (2003) IMF values from the Skel¬ 
ton et al. (2014) catalogue by a factor of 1.8. We divide 
the sample of photo-z galaxies at z = 0 — 4 into two 
subsamples of star-forming galaxies (SFGs) and quies¬ 
cent galaxies (QGs) by the rest-frame UVJ color criteria 
of Muzzin et al. (2013). Because the UVJ color criteria 
are not tested for z > 4 sources, we do not apply these 
color criteria to the photo-z galaxies at z > 4. Muzzin 
et al. (2013) find that the QG fraction is small, 10%, 
at z ~ 3.5, and it is likely that a QG fraction at the 
early cosmic epoch of z > 4 is negligibly small, perhaps 
< 10%. We thus regard all of the z > 4 photo-z galax¬ 
ies as SFGs. The total numbers of SFGs and QGs are 
165,517 and 10,631, respectively. The iJieo magnitude 
at the 50% completeness is ~ 26.5 mag for the photo- 
z galaxies in deep CANDELS fields. The details of the 
completeness estimates and values are presented in Skel¬ 
ton et al. (2014). 

2.2. Sample of LBGs at z = 4 — 10 
in CANDELS, HUDF09/12, and HFF 

The second sample consists of 10,454 LBGs at z = 
4 — 10 made by Y. Harikane et al. (in preparation) in the 
CANDELS, the Hubble Ultra Deep Field 09-bl2 (HUDF 
09-1-12; Beckwith et al. 2006; Bouwens et al. 2011; Illing¬ 
worth et al. 2013; Ellis et al. 2013) fields®, and the par¬ 
allel fields of Abell 2744 and MACS0416 in the Hubble 
Frontier Fields (e.g., Coe et al. 2014; Atek et al. 2015; 

® http://3dhst.research.yale.edu/Home.html 

® http://archive.stsci.edu/prepds/xdf/ 


Oesch et al. 2014; Ishigaki et al. 2014). The numbers 
of our LBGs are summarized in Table 3. These LBGs 
are selected with the color criteria, similar to those of 
Bouwens et al. (2014b). We perform source detections by 
SExtractor (Bertin & Arnouts 1996) in coadded images 
constructed from bands of Yo 98 Uio 5 '/i 25 -ffi 60 , J 125 H 1 G 0 , 
and Higo for the z ^ 4 — 7, 8, and 10 LBGs, respectively. 
The JHi 4 o band is included in the coadded image for 
the z - 7 - 8 LBGs in the HUDF09-hl2 field. The flux 
measurements are carried out in Kron (1980)-type aper¬ 
tures with a Kron parameter of 1.6 whose diameter is 
determined in the iLieo band. In two-color diagrams, we 
select objects with a Lyman break, no extreme-red stel¬ 
lar continuum, and no detection in passbands bluer than 
the spectral drop. See Y. Harikane et al. (in prepara¬ 
tion) for more details of the source detections and LBG 
selections. 

The Higo magnitudes at the 50% completeness is ~ 28 
mag for the LBGs in deep CANDELS fields (Bouwens 
et al. 2014b). The details of the completeness estimates 
and values are presented in Y. Harikane et al. (in prepa¬ 
ration) . 

Several previous studies on galaxy size have included 
a galaxy at z ^ 12 selected in photo-z technique (Ellis 
et al. 2013). In this study, we do not use the galaxy at 
z ~ 12 because the redshift of the source is under debate 
(e.g., Ellis et al. 2013; Brammer et al. 2013; Bouwens 
et al. 2013; Capak et al. 2013; Pirzkal et al. 2013). 

2.3. Stellar Masses of Photo-z Galaxies and LBGs 

Some analyses and discussions in this work require M* 
of the photo-z galaxies and the LBGs. For the photo-z 
galaxies, we take values from Skelton et al. (2014). 
For the LBGs, we derive stellar masses, adopting an em¬ 
pirical relation between UV magnitude Muv and M*. 
First, we calculate Muv from the total magnitudes in 
the LBG detection images (Section 3), assuming that the 
typical redshifts are (z) ~ 3.8,4.9, 5.9, 6.8, 7.9, and 10.4. 
The stellar masses are obtained by converting their Muv 
through the empirical Gonzalez et al.’s relation (see also, 
the updated result of Gonzalez et al. 2014), 
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Fig. 2.— Examples of Sersic profile fitting results. The right, middle, and left panels indicate the results for the SFGs in the rest-frame 
optical and UV stellar continuum emission, and the LBGs, respectively. Each of 4-panel sets presents the original, the best-fit model, the 
residual, and the mask images, from left to right (see Section 3 for details) for one object. Each row, from top to bottom, denotes example 
galaxies from 2 : 0 to 10. The SFGs at 2 : = 1 — 2 and 2 : = 2 — 3 are the same objects between the left and right panel sets that exhibit 
the images at the rest-frame optical and UV wavelengths, respectively. The black tick in the panel for the SFG at 2 : 0 — 1 indicates the 

size of 1". North is up and east is to the left. 


TABLE 4 

Catalog of Photo-2: Galaxies with our Size Measurements 


Catalog ID 

muv 

RUV . 
e.maior 

nuv 

<?uv 

mopt 

pOpt 

e,major 

^Opt 

90pt 

flag 


[mag] 

[arcsecj 



[mag] 

[arcsec] 




(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

gds_29 





23.70 ±0.021 

0.489 ± 0.016 

0.94 ±0.04 

0.92 ± 0.02 

0 

gds_32 





24.34 ± 0.042 

0.164 ±0.013 

4.15 ±0.35 

0.77 ± 0.03 

0 

gds_59 





25.31 ±0.029 

0.183 ± 0.008 

0.65 ±0.09 

0.57 ±0.03 

0 

gds_86 





25.58 ±0.031 

0.174 ±0.008 

0.36 ±0.09 

0.73 ± 0.04 

0 

gds_122 





25.70 ±0.051 

0.201 ± 0.017 

0.66 ±0.14 

0.87 ± 0.06 

0 


Note. — A catalog of the photo -2 galaxies with S/N > 15 and reliable outputs of GALFIT fitting. Five example objects 
are shown here. Columns: (1) Catalog ID. The alphabetical characters represent the HST fields (“gds”; GOODS-South, 
“gdn”; GOODSN-North, “uds”; UDS, “aeg”; AEGIS, “cos”; COSMOS). The numeric characters correspond to the ID 
number in the 3D-HST catalog (Skelton et al. 2014). (2) and (6) Total magnitude. (3) and (7) Effective radius along 
the major axis in arcseconds. (4) and (8) Sersic index. (5) and (9) Axis ratio. (10) Flag for the reliability of GALFIT 
fitting. The values of 0 and 1 indicates reliable and unreliable measurements, respectively. (2)-(5) Measurements at 
— 1500 — 3000 A. (7)-(9) Measurements at — 4500 — 8000 A. All measurement uncertainties are the half-width 

of the 68%-confidence interval. 

(The complete table is available in a machine-readable form in the online journal.) 


logM*— —39.6 + 1.7 X logL i5oo, , , 
= -4.50-0.68 X Muv, 

where ii 5 oo is the luminosity at the rest-frame 1500 A. 
This empirical relation is derived under the assumptions 
similar to ours (the Salpeter IMF and no nebular emis¬ 
sion lines included in SED). 

To test whether this empirical relation (eq. 1) of TuV" 
M, is reliable and consistent with the M* estimates of 


the photo-z galaxy sample, we compare this empirical 
relation with Muy-M, relations derived from the photo- 
z galaxies. 

We estimate Muv from the absolute UV magnitudes at 
a wavelength of 2800 A from the photo-z catalog, assum¬ 
ing the majority of star-forming galaxies have a flat UV 
spectrum of fv =const. We present Muy-M, relations 
of the photo-z galaxies in Figure 1. The UV magnitude 
correlates well with M*, suggesting the existence of the 
“star-formation main sequence” (e.g., Daddi et al. 2007; 
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TABLE 5 

Catalog of LBGs with our Size Measurements 


Catalog ID 

muv 

RUV . 
e,major 

Quv 

flag 


[mag] 

[arcsec] 



(1) 

(2) 

(3) 

(4) 

(5) 

z4_gdsd_10028 

26.42 ±0.05 

0.344 ± 0.029 

0.16 ± 0.03 

0 

z4_gdsd_10045 

25.23 ±0.02 

0.241 ± 0.008 

0.57 ±0.02 

0 

z4_gdsd_10054 

26.58 ±0.03 

0.100 ±0.009 

0.47 ± 0.08 

0 

z4_gdsd_10153 

27.60 ±0.07 

0.102 ±0.024 

0.43 ± 0.21 

0 

z4_gdsd_10202 

26.33 ±0.04 

0.196 ±0.012 

0.38 ± 0.05 

0 


Note. — A catalog of the LBGs with S/N > 15 and reliable outputs 
of GALFIT fitting. Five example objects are shown here. Columns: (1) 
Catalog ID in Y. Harikane et al. (in preparation) (2) Total magnitude. 
(3) Effective radius along the major axis in arcseconds. (4) Axis ratio. 
(5) Flag for the reliability of GALFIT fitting. The values of 0 and 1 
indicates reliable and unreliable measurements, respectively. The GALFIT 
fitting is performed in the coadded HST images (see Section 3). Note 
that Sersic indices are not listed due to fixed n values in the GALFIT 
fitting. All measurement uncertainties are the half-width of the 68%- 
confidence interval. 

(The complete table is available in a machine-readable form in the online 
journal.) 
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Fig. 3.— Comparison between the effective radii of re and r-e^vdw 
measured by this study and van der Wei et al. (2014), respectively, 
for objects with S/N > 15 (black dots) and < 15 (gray crosses 
at Hiqo 26.5 mag). The right panel shows histograms for the 
number of the galaxies. The black and gray histograms denote 
objects with S/N > 15 and < 15, respectively. The number of 
objects with S/N < 15 is multiplied by a factor of 100 for clarity. 


Lee et al. 2012; Whitaker et al. 2012; Steinhardt et al. 
2014). 

Figure 1 presents that the slopes of the relations appear 
to be flatter at a bright range of Muv ^ ~22 than at a 
faint Muv range. Similar flat slopes are reported by a 
large survey area of the CANDELS fields (Stark et al. 
2009; Lee et al. 2011; Salmon et al. 2015). Because our 
LBGs used in this analysis have magnitudes of Muv > 
—22, we fit logM, = a + 6Muv to the Muv-M* relation 
at Muv > —22, where a and b are free parameters. The 
best-fit functions for the photo-z galaxies are 


log [Mq] = 


1.46-0.43 X Muv(2 
0.82 — 0.45 X Muy[z 
1.12-0.43 X Muy{z 
—0.22 — 0.49 X M\jy(^z 
-2.10-0.58 X Muv(^ 
-2.45- 0.59 X Muv(z 


0 

1 

2 

3 

4 

5 


1 ), 

2 ), 

3 ) , 

4 ) , 

5 ) , 

6 ) . 


(2) 


If we assume that the magnitudes of 1700 A to 1500 A 
are the same for typical LBGs with f^, =const, the slopes 
b of —0.58 ± 0.02 and —0.59 ± 0.03 at z ^ 4 — 6 roughly 
agree with that of eq. (1) (i.e., b = —0.68±0.08). We thus 
conclude that the empirical relation (eq. 1) is reliable and 
consistent with the M* estimates of the photo-z galaxy 
sample. Moreover, no strong evolution in the Muv-M* 
relation is found at z > 4 in eq. (2) and Figure 1. We 
use eq. (1) to estimate M* of our z > 4 LBGs. 


3. SIZE MEASUREMENT 


In this section, we describe methods to measure galaxy 
sizes by using the high spatial resolution images oiHST. 
To minimize the effect of morphological A-correction, we 
use images of four bands, Vgoe and /814 on ACS^, and 
J 125 and Hieo on WEC3/IR. We select one of these bands 
whose entire passband is covered by the wavelength range 
of = 1500 - 3000 A or A°f = 4500 - 8000 A of 
each object. If two or more filter passbands meet this 
criterion, we chose a band that observes the shortest 
wavelength. Prior to the size measurements, we extract 
18" X 18" cutout images from the data 

at the position of each photo-z galaxy and LBG. The size 
of cutout images is sufficiently large to investigate entire 
galaxy structures even for extended objects at z ^ 0 — 1. 
We use coadded images of I 098 L 105 T 125 L 7160 , ^ 125 ^^ 160 , 
and iLieo constructed in Section 2.2 for the z ^ 4 — 7, 8, 
and 10 LBGs, respectively. The limiting magnitudes of 
these coadded images are summarized in Table 1. 

We measure the galaxy size basically in the same man¬ 
ner as previous studies for high-z LBGs (e.g., Ono et al. 
2013) based on the two-dimensional (2D) surface bright¬ 
ness (SB) profile fitting with the GALFIT software (Peng 
et al. 2002, 2010). We fit a single Sersic profile (Sersic 
1963, 1968) to the 2D SB distribution of each galaxy to 
obtain the half-light radius along the semi-major axis, 
■Ce,major- The major IS Converted to the “circularized” 

radius, Tg, through re = a^/bja = re.major-/?! where a, 
6, and q are the major, minor axes, and axis ratio, re¬ 
spectively. Several authors studying z 0 — 3 galaxies 
claim that Ce,major should be used, because rg,major does 
not depend strongly on the galaxy inclination (e.g., van 
der Wei et al. 2014). However, the circularized radius 
Te has widely been used in size measurements for faint 
and small high-z sources (e.g., Mosleh et al. 2012; Ono 
et al. 2013; Holwerda et al. 2014). We here use the circu¬ 
larized radius re in order to perform self-consistent size 
measurements and fair comparisons from z ~ 0 to ~ 10. 

We create sigma and mask images for estimating the 
fitting weight of individual pixels and masking neighbor¬ 
ing objects of the main galaxy components, respectively. 


^ We make use of zsso for GOODS-North that has not been 
taken with the /814 band. 
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Fig. 4.— Difference between and for the SFGs at 1.2 < 
z < 2.1 (gray dots) as a function of stellar mass. The gray circles 
represent the median values of in different stellar 

mass bins. The right panel shows histograms for the number of 
SFGs. 


The sigma images are generated from the drizzle weight 
maps produced by the ffST data reduction (Koekemoer 
et al. 2003). We also include the Poisson noise from 
the galaxy light to the sigma image (e.g., Hathi et al. 
2009; van der Wei et al. 2012). The mask images 
are constructed from segmentation maps produced by 
SExtractor. We identify neighboring objects with the 
SExtractor detection parameters of DETECT_MINAREA= 
5 pixel, DETECT_THRESH= 2 ct, DETECT_NTHRESH= 16, and 
DEBLEND J1INC0NT= 0.0001. 

We input initial parameters taken from the 3D- 
HST+CANDELS photometric catalogue (Skelton et al. 
2014) for the photo-z galaxies. Specifically, the total 
magnitude m, axis ratio q, position angle P.A., and half 
light radius re of each galaxy are initial parameters that 
are written in the GALFIT configuration file. The Sersic 
index n is set to n = 1.5 as an initial value for the photo- 
z galaxies, while initial n does not affect strongly fitting 
results (Yuma et al. 2011, 2012). In fact, we change the 
initial parameters of Sersic index to n = 1 and 3, but still 
obtain similar best-fit n values even with these different 
initial parameters. For the LEG sample, the initial pa¬ 
rameters are taken from the results of SExtractor pho¬ 
tometry (Y. Harikane in preparation). The Sersic index 
for LBGs is fixed to 1.5 for reliable fitting for faint and 
small high-z sources. This fixed Sersic index is justified 
by the evolution of n in SFGs, as demonstrated in Sec¬ 
tion 5.1. To obtain re, n, and q, we allow the parameters 
to vary in the ranges, Am < 3 mag, 0.3 < Vg < 400 
pixels, 0.2 < n < 8, 0.0001 < g < 1, Ax < 2 pixel, and 
Ay < 2 pixel, which are quite similar to those of van 
der Wei et al. (2012). We discard objects whose one or 
more htting parameters reach the limit of the parameter 
ranges (e.g., Ve = 400). The PSF models of the HST 
images are provided from the 3D-HST project (Skelton 
et al. 2014). 

We have analyzed the photo-z galaxies and LBGs 
shown in Sections 2.1 and 2.2. As we discuss below, sizes 
of faint galaxies are poorly determined. We thus choose 


photo-z galaxies and LBGs whose sources have a signal- 
to-noise ratio (S/N) greater than 15. This S/Athreshold 
is determined by Monte Carlo simulations for faint and 
small high-z sources (e.g., van der Wei et al. 2012; Ono 
et al. 2013). Tables 2 and 3 summarize the number of 
photo-z galaxies and LBGs, respectively, that are ana¬ 
lyzed in our study. The object numbers in our size anal¬ 
ysis are 142,273 (9,767) in Vgoe, 136,493 (10,118) in / 814 , 
139,308 (10,845) in J 125 , and 147,204 (11,297) in Hiqq, 
for the SFGs (QGs) of the photo-z sample, and 7,233 
for the LBGs. The total numbers of SFGs (QGs) that 
are well fit in the optical and UV stellar continuum emis¬ 
sion are 59,115 (4,234) at z ~ 0 - 3 and 30,765 (799) 
at z ^ 1 — 6, respectively, while sizes of 4,993 LBGs are 
securely measured. Tables 4 and 5 show the size mea¬ 
surements given by our structural analysis for the photo- 
z galaxies and the LBGs, respectively. Figure 2 presents 
example images of the fitting results, demonstrating that 
our size measurements are well performed. 

Note that clumpy structures are masked in the fitting, 
as indicated in the mask panels of Figure 2. This mask¬ 
ing procedure is included in our analyses, because a sin¬ 
gle Sersic profile fitting is not reliable for galaxies with 
the clumpy structures. Moreover, the number of well-fit 
galaxies decreases, if no masking is applied. Neverthe¬ 
less, we examine whether the masking procedures change 
our conclusions, and find that the Vg measurements are 
statistically comparable in galaxies with and without 
masking. The fraction of galaxies with the clumpy struc¬ 
tures ranges from ^ 30% at z ~ 1 to ^ 50% at z ~ 2. 
This study only addresses galaxies’ major stellar compo¬ 
nents. The detailed analyses and the results of clumpy 
stellar sub-components are presented in the paper 11. 

van der Wei et al. (2012, 2014) obtain their re val¬ 
ues in the J 125 , JH 14 Q, and iAieo bands for all the 3D- 
HST-I-GANDLES galaxies using the GALAPAGOS software 
(Barden et al. 2012) which is a wrapper of SExtractor 
and GALFIT for morphological analyses. Several mor¬ 
phological studies have utilized GALAPAGOS allowing for 
the simultaneous determination of both the structural 
parameters and the background flux level for multi¬ 
objects. In Figure 3, we compare rg measurements of 
ours with those of van der Wei et al. (2014) estimated 
with GALAPAGOS. We find that our rg values are in good 
agreement with those obtained by van der Wei et al. 
(2014). We also find that faint galaxies with S/N < 15 
are significantly scattered in Figure 3. This confirms that 
the threshold of S/N > 15 is important for secure size 
measurements. 

4. iL-CORRECTION, STATISTIGAL CHOIGE, AND 
SELEGTION BIAS 

4.1. Effect of Morphological K-Correction 

We investigate the effects of morphological K- 
correction in our size measurements, comparing our rg 
at different wavelengths. Because the HST imaging data 
covers up to iAieo band, we can study the rest-frame 
UV morphology for galaxies at z > 3. Understanding 
the effects of morphological A-correction is considerably 
important to evaluate the size evolution of star-forming 
galaxies over a wide redshift range of 0 < z < 10. The 
sizes in the rest-frame UV and optical stellar continuum 
emission, and tracing different stellar pop- 
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TABLE 6 

Summary of the Best-Fit Size Growth Rates 


Data points 

(1) 

Sample 

(2) 

LuvILl^^ 

(3) 

[kpc] 

(4) 

Hz 

(5) 

Bh 

[kpc] 

(6) 

Hh 

(7) 

Median 

All 

1-10 

0.3-1 

0.12-0.3 

4.78 ± 0.68 
5.45 ±0.31 
4.44 ±0.19 

-0.84 ±0.11 
-1.10 ±0.06 
-1.22 ±0.05 

3.80 ± 0.40 
4.33 ±0.17 
3.46 ±0.13 

-0.62 ± 0.08 
-0.86 ± 0.04 
-0.97 ±0.05 


w/o rPP* 

1-10 

0.3-1 

0.12-0.3 

4.05 ±0.59 
5.21 ±0.28 
3.54 ±0.58 

-0.78 ±0.08 
-1.15 ±0.07 
-1.11 ±0.11 

3.09 ± 0.36 
3.54 ±0.29 
2.45 ± 0.32 

-0.56 ± 0.06 
-0.80 ±0.05 
-0.78 ±0.08 

Average 

All 

1-10 

0.3-1 

0.12-0.3 

5.80 ±0.65 
5.85 ±0.33 
5.52 ±0.43 

-0.79 ±0.10 
-0.95 ±0.07 
-1.17 ±0.07 

4.91 ±0.42 
4.83 ±0.20 
4.29 ±0.27 

-0.61 ± 0.07 
-0.74 ± 0.04 
-0.87 ±0.05 


w/o rp^* 

1-10 

0.3-1 

0.12-0.3 

11.3 ±4.44 
10.9 ±2.94 
6.82 ± 2.25 

-1.22 ±0.25 
-1.36 ±0.18 
-1.31 ±0.20 

7.48 ± 2.37 
6.90 ± 1.50 
4.37 ± 1.18 

-0.85 ± 0.18 
-0.95 ± 0.13 
-0.91 ± 0.14 

Mode 

All 

1-10 

0.3-1 

0.12-0.3 

4.00 ± 0.49 

4.45 ± 0.89 
3.28 ±0.18 

-0.78 ±0.08 

-1.26 ±0.17 
-1.23 ±0.07 

3.07 ±0.30 

2.97 ±0.45 
2.56 ±0.11 

-0.55 ± 0.57 

-0.89 ± 0.12 
-1.00 ±0.05 


w/o rp^* 

1-10 

0.3-1 

10.9 ±3.94 
3.00 ±0.19 

-1.14 ±0.25 
-1.01 ±0.05 

7.45 ±2.14 
2.15 ±0.10 

-0.80 ±0.17 
-0.71 ± 0.03 


0.12-0.3^^ 

Note. — Columns: (1) Statistics of re- (2) Sample used in the fits for the size evolution. “All” denotes 
to use all samples in a Luv bin. “w/o represents to exclude the data points of (3) Bins 

of Luv iri units of L*_^. (4) of - 6^(1 + . (5) of ^^{1 + . (6) Bfj of B}ih{z )^^, where 

h(z) = H{z)/Ho = vTuITT^F+nl. (7) I3h of BHhizfH, 

^ The minimization is not converged. 

ulation would yield a large difference in Vg. Here we 
make a comparison between and of the SFGs 
at 1.2 ^ z < 2.1 where the both radii can be measured 
with the HST data. 

Figure 4 shows the differences between and of 
the SFGs as a function of stellar mass. Although we find 
a large scatter, the median values of (r^'^ — r®P*)/r^P* 
are less than 20 % in all stellar mass bins. This indicates 
that the differences of statistical rg measurements are 
small for star-forming galaxies with logM* =9—11 Mq 
at z 1 — 2. 

Similarly, van der Wei et al. (2014) have found that 
?’e,major is typically smaller in redder bands for SFGs at 
z ^ 0 — 2 (see also, e.g., Szomoru et al. 2011; Wuyts 
et al. 2012). This trend is more significant in more mas¬ 
sive SFGs. The smaller size in redder bands would be 
interpreted as heavier dust attenuation in the galactic 
central regions in bluer bands (e.g., Kelvin et al. 2012) 
and/or the inside-out disk formation (e.g., Brooks et al. 

2009; Bezanson et al. 2009; Naab et al. 2009; Nelson et al. 

2012; Patel et al. 2013). We confirm the wavelength de¬ 
pendence even in our rg in the most massive M* bin 
as shown in Figure 4. van der Wei et al. (2014) have 
parametrized the wavelength dependence of re^major as a 
function of redshift and stellar mass. Following the for¬ 
mula, the size difference fraction (f^^ —r°P*)/rpP* is cal¬ 
culated to be c." 30% for z ~ 2 galaxies with logM* = 11 
Mq. 

Note that the difference of stellar population becomes 
smaller at z > 2 than z ^ 1 — 2, because the short cos¬ 
mic age of z > 2 provides a smaller stellar-age difference 
and a less metal enrichment than that of z ~ 1 — 2. 

This agreement of and r^P* suggests that the sta¬ 
tistical values represent the typical sizes of stellar- 
component distribution for star-forming galaxies of SFGs 
and LBGs at z > 3 with a small systematic uncertainty 


of <30 %®. 

We examine the effect of morphological A-correction 
in more detail by investigating evolutionary trends of Vg 
and size-relevant quantities in the rest-frame optical and 
UV emission for the photo-z galaxies at z ^ 0 — 6. Figure 
5 presents redshift evolution of rg, n, and star-formation 
rate surface density (SFR SD), Esfr- The SFR SD is 
derived in the effective radius, and calculated by 

SsFR [Mq yr“^ kpc“^] = , (3) 

7rr/ 

where a factor of 1/2 corrects for the SFR value which 
is derived from the total magnitudes. For the photo-z 
galaxies, we use SFRs taken from the catalog of Skelton 
et al. (2014). For LBGs, we compute SFRs from Lyv 
using the relation of Kennicutt (1998a), 

SFR [Mq yr-p = 1.4 x lO’^^L^ [ergs“^ Hz"p. (4) 

van der Wei et al. (2014) have already examined the 
Tg evolution at 5000 A in the rest-frame for galaxies at 
0 ^ ^ 3 in the 3D-HST-I-CANDELS sample. In our 

study, we extend this analysis of0<z<3toz>4, 
using the photo-z galaxies and the LBGs. 

In Figure 5, median values of these quantities are in 
good agreement between the measurements in the rest- 
frame optical and UV emission of the SFGs at 2 < z < 3. 
Additionally, the evolutionary tracks at z < 3 smoothly 
connect with those at z > 3. We also find no strong 

® In Figure 4, we find that the scatters of 
are comparably large in high and low-mass galaxies. Because the 
scatters originated from statistical errors should be smaller in the 
high-mass galaxies than the low-mass galaxies, the scatters of the 
high-mass galaxies are probably not dominated by statistical errors 
but intrinsic re differences. 
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Fig. 5.— Redshift evolution of size-relevant quantities for the photo-z galaxies in different stellar mass bins (logM* = 9 — 9.5, 9.5 — 
10,10 — 10.5, and 10.5 — 11 Mq from left to right). From top to bottom, the panels show the effective radius Sersic index n, and SFR 
SD EgpQ. The blue and red circles indicate median values of and respectively, for the SFGs. The open squares in the rightmost 

panels present median values of r®*’* for the QGs with log M, = 10.5 — 11 Mq. The error bars denote the 16th and 84th percentiles of the 
data point distribution. In the top panels, the best-fit re curve in the bin of logM* = 9.5 — 10 Mq is plotted for reference with the solid 
and dashed lines. The best-fit and /3ff values are —0.72 ± 0.04 and —0.60 ± 0.02, respectively. The dashed gray lines in the bottom 
panels represent the SgpG evolution calculated with SFR= lOMQ/yr and the best-fit re curve. The horizontal lines in the second-top 
panels denote a weighted mean of (n) = 1.4. In the second-right and right-most panels of the Sersic index plots, the open diamonds denote 
measurements of Sersic index for SFGs in Morishita et ah (2014). The Sersic index estimates of Morishita et ah (2014) are comparable to 
ours. In the panels of low mass log M* = 9 — 9.5 Mq, the quantities for high redshifts are not plotted, due to their poor statistics. 


dependence of these evolutionary trends on the stellar 
mass. These agreements confirm a small effect of mor¬ 
phological iGcorrection in the median re values. 

4.2. Statistical Difference and Selection Bias 

We examine redshift evolution of median, average, and 
modal re of our galaxies to evaluate statistical differences 
and selection biases. We define four Lyy bins for these 
analyses. The Luv-bins are 1 — 10, 0.3 — 1, 0.12 — 0.3, 
and 0.048 — 0.12 Luv/-^z= 3 j where is the charac¬ 
teristic UV luminosity of LBGs at z ^ 3 (Myv = — 21, 
Steidel et al. 1999)®. To investigate the re distribution 
shape, in Figure 6 we plot the re distribution of SFGs 
and LBGs at z ^ 1 — 6 in the bin of Luv = 0.3 — 1 
that has good re measurement accuracies whose typical 
reduced values are the smallest among the Tyv bins of 

® The Luv-bins are the same as in previous studies (e.g., Oesch 
et al. 2010). The LBGs in the faintest Lyv bin are used only for 
the stacking analysis (Section 5.2). 


the re measurements. We fit the re with the log-normal 
distribution. 


Pi^e) 


1 

- i= exp 

^e^ln re V 


\T?{re/re) ' 

2(7? 

In re -I 


( 5 ) 


where ry and uinre are the peak of re and the standard 
deviation of In re, respectively. We fit the log-normal 
functions to the re-distribution data with two free pa¬ 
rameters of ry and uinrej and present the best-ht log¬ 
normal functions in Figure 6 for the data of good statis¬ 
tics, the SFGs at z ~ 1 — 3 and the LBGs at z ~ 4 — 6 
in the Lyv = 0.3 — l.OLyy bin. The re distributions 
of the high-z star-forming galaxies are well represented 
by the log-normal distribution. The reduced values 
are 0.006, 0.003, 0.004, 0.005, and 0.011 for the SFGs at 
z = 1 — 2 and 2 — 3, and the LBGs at z ^ 4, 5, and 6, 
respectively. Figure 7 is the same as Figure 6, but for all 
of our galaxies. Figure 7 indicates that the rg distribu- 
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TABLE 7 

Summary of the LBG Size Growth Rates from the Previous Studies 


Reference 

Number 

Redshift Range 

I3z of (1 -b z)'3- 

Statistics 

Size Measurements 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

Bouwens et al. (2004) 

■ ■ ■ (2929) 

2-6 

-1.05 ±0.21 

Average 

SExtractor 

Ferguson et al. (2004) 

■ • ■ (773) 

2-5 

-1.5 

Average 

SExtractor 

Ravindranath et al. (2006) 

1333 (4694) 

3-5 



GALFIT 

Hath! et al. (2008a) 

61 (61) 

3-6 

-1.5 

Average 

SExtractor 

Conselice S>c Arnold (2009) 

583 (583) 

4-6 



SExtractor 

Oesch et al. (2010) 

21 (21) 

7-8 

-1.12 ±0.17 

Average 

SExtractor, GALFIT 

Grazian et al. (2012) 

(153) 

7 



SExtractor 

Mosleh et al. (2012) 

■■■ (218) 

4-7 

-1.20±0.11 

Median 

GALFIT 

Huang et al. (2013) 

1012 (1356) 

4-5 

r'../ —1 

Mode 

SExtractor, GALFIT 

Ono et al. (2013) 

15 (81) 

7-10 

1 

CO 

o 

1 + 

o c 

Average 

GALFIT 

Curtis-Lake et al. (2014) 

1318 (3738) 

4-9 

-0.31 ±0.26 

Mode 

SExtractor 

Holwerda et al. (2014) 

8 (8) 

9-10 

-1.0 ±0.1 

Average 

GALFIT 

Kawamata et al. (2014)^ 

39 (39) 

6-8 

-1.24 ±0.1 

Average 

glafic 

This work 

4993 (10454) 

4-10 

-1.10 ±0.06 

Median 

GALFIT 




-0.95 ±0.07 

Average 

GALFIT 




-1.26±0.17 

Mode 

GALFIT 

inch Photo- 2 : SFGs 

89880 (312722)'= 

0-6 





Note. — Columns: (1) Reference. (2) Number of galaxies whose size is measured in the reference. The values in parentheses 
are the number of galaxies in parent sample. (3) Redshift range for size measurements of LBGs. (4) Best-fit pz of (1 + 2 )^ 
for a bright (Luv ~ 0.3 — 1 L*_^) galaxy sample. (5) Statistics for deriving a representative rg at a redshift. “Mode ” 
corresponds to the peak of size distribution derived by the fitting with a log-normal function (Equation 5). (6) Method or 
software to measure galaxy sizes. 

^ Sample galaxies are selected in a field of galaxy cluster. This study corrects for the gravitational lensing effects of magnifi¬ 
cation and shear with their mass model. 

^ The value is a total number of SFGs whose sizes are well measured in and See Table 2. 


tions are well fitted by the log-normal functions in the 
wide ranges of redshift, z ~ 0 — 6, and the UV luminosity, 
~ 0.12 — lOLyy. Note that log-normal functions cannot 
be fitted to the data of the z>7 galaxies and some low-z 
galaxies in Figure 7, due to the small statistics. More¬ 
over, the fitting result of z ^ 0—1 is only obtained for the 
distribution in the luminosity bin of 0.12 — O.lLJjy 
because of the poor statistics of the other-luminosity bin 
data. 

Because the re-distributions follow the log-normal 
functions, the average, median, and modal values of In re 
should be the same. However, in the previous studies, 
the size evolution is discussed with the average, median, 
and modal values of in the linear space (e.g., Hol- 
werda et al. 2014; Ono et al. 2013; Grazian et al. 2012; 
Oesch et al. 2010; Bouwens et al. 2004). Here we obtain 
re measurements with different statistics choices in the 
linear space, following the previous studies, and evaluate 
the differences of the size evolution results. We derive 
size growth rates based on average, median, and modal 
re in a bin of 0.3 — 1 Lijy /estimating the modal 
by fitting size distributions with a log-normal function. 
In Figure 8, we compare our measurements with those 
of the previous studies that apply the different statistics. 

We confirm that our results are consistent with those 
of the previous studies. Moreover, Figure 8 indicates that 
galaxy sizes decrease from z ^ 0 to ^ 6 in any statistical 
choices of average, median, and mode. 

We fit rg = Bz{l + z)^* for the average, median, 
and modal rg values given by our and previous stud¬ 
ies, where and j3z are free parameters. The fitting 
is performed for the combination of and as 

Because there are only three LBGs at z ~ 10, the weighted 
average rg is only derived for our z ~ 10 LBGs. Note that the 
z ~ 10 data is presented in Figure 8, but that the data is not used 
to derive the size evolution function below. 


well as for only. Table 6 summarizes the best-fit 
Bz and Pz values. Table 7 is a summary of the samples 
and f5z values from our and previous studies for LBGs 
z > 4. Our average, median, and modal rg values scale 
as oc (1 -I- z)~^ indicating that, again, the choices 
of statistics in rg measurements give no significant im¬ 
pacts on size growth rates. This conclusion is consistent 
with the result that shows no significant evolution 

as discussed in Section 6.1.1. 

Most previous studies have employed average values 
for representative rg. However, Figure 7 indicates that 
the median measurements trace the typical galaxy sizes 
parametrized by r7 better than the average values. Be¬ 
cause the small samples of z > 7 galaxies do not allow 
us to estimate modal rg values, we use median values for 
our main analyses, unless otherwise specihed. 

Figure 8 compares the rg values of SFGs and LBGs at 
z ^ 4 — 6 with a bright UV luminosity. In any statistics 
choices, we find that the rg values of SFGs and LBGs are 
comparable within the scatters of < 30%. These results 
indicate that star-forming galaxies selected by photo-z 
and dropout techniques statistically give the similar rg 
values, and that the bias from the different selection tech¬ 
niques is as small as < 30% in the rg determination. 

5. RESULTS 
5.1. Sersic Index 

A Sersic index represents the SB profiles of galaxies. 
A high n means a cuspier SB distribution, indicating the 
existence of a central bulge. On the other hand, a lower 
n suggests a disk-like light profile with a flatter SB dis¬ 
tribution at the central galactic region. The Sersic index 
depends on observed wavebands and stellar populations 
(e.g., color), which have been revealed by detailed struc¬ 
tural analyses with multiple passbands for local galaxies 
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TABLE 8 

Size-Luminosity Relation at z = 0 — 8 


Muv 

,,Opt 

' e 

Muv 

^uv 

Muv 

^uv 

[mag] 

[kpc] 

[mag] 

[kpc] 

[mag] 

[kpc] 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

z = 0 

- 1SFGs 

2 = 1 

- 2 SFGs 


4 LBGs 

-21.0 


-21.0 


-23.0 

1.536lJ;°" 

-20.0 

3.380l?'^^f 

-20.0 

1-79412;"® 

-22.0 

1 ooq+l-lOO 
l.ZZO_Q QQQ 

-19.0 

Z.OU4_;L.045 

-19.0 

1 r)q9+l-464 

i.zoz_Q 034 

-21.0 

1.0581®;®" 

-18.0 

1 Qq7’+i-234 

1.00f_Q gJQ 

-18.0 

0.940lJ°®« 

-20.0 

0.7331®;®^ 

-17.0 

1 O4Q+0.983 
1.C)4U_q 5YQ 

-17.0 

0.723lJ;°i® 

-19.0 

0.5891®;"° 

-16.0 

l-074l“;Si 

-16.0 

U.DUD_o.404 

-18.0 

0.509l®;«« 

-15.0 

0.854/“;8|° 

-15.0 

n 461 +0-398 
G.^01 _o.309 

-17.0 

0.4381®;®®® 

2: = 1 

- 2 SFGs 

z = 2 

- 3 SFGs 

2 

5 LBGs 

-22.0 

q nqK+0.237 

O.UOO_Q 44g 

-22.0 

1 z19S+^-002 

-*-•^^^-0.632 

-22.0 

i.025lj;§^^ 

-21.0 

1.958l?;“4 

-21.0 

-‘-•+^'^-0.727 

-21.0 

0.788l®;i® 

-20.0 

1.982lJ:«4 

-20.0 

1 n'7£j+l-177 
1-076Ios4s 

-20.0 

0.595l®;i® 

-19.0 

l.c>OU_Q 5gg 

-19.0 

n 81 q+O-787 
U.OiO_o.382 

-19.0 

0.519l®;«6 

-18.0 


-18.0 

0.685l®;J« 

-18.0 

U.OUD_o.284 

-17.0 

1.0281^1? 

-17.0 

0.5091°;®^® 

-17.0 

n ‘4fyG+0 095 

U.OOU_o.o76 

-16.0 

0-937+J:!^/^ 

z = 3 

- 4 SFGs 

2 

6 LBGs 

z = 2 

- 3 SFGs 

-22.0 

1 470+1-687 
1'4<o_o.705 

-22.0 

i.uoo_o.696 

-21.0 

Q o'70”I”':5-^>19 / 

z.o<o_Q 908 

-21.0 

i-054lj;i®^ 

-21.0 

0.635l®;"I 

-20.0 

1 9c;q+0.488 

i.zoo_o.o98 

-20.0 

0-778l°;I®i 

-20.0 

0.5651®;™° 

-19.0 

1.UOD_o.7io 

-19.0 

0.6201®;®®® 

-19.0 

0.584l®;"4 

-18.0 

1.0ZO_o.9g7 

-18.0 

u.o/z_0 269 

-18.0 

n “^71 +0-222 

U.O/ -L_o.198 



2 = 4 

- 5 SFGs 

2 

7 LBGs 



-22.0 

i-08i1^;2®s 

-21.0 

n 7q7+0.320 
-0.421 



-21.0 

0.892l°;^e? 

-20.0 

u.^oy_o 268 



-20.0 

0.7081®;™® 

-19.0 

U.4iO^_o 207 



-19.0 

0 444+0-362 

2 

8 LBGs 



z = 3 

- 6 SFGs 

-21.0 

0.419lJ;l®i 



-22.0 

0 975+'^-''^'' 
u.y /O_o 425 

-20.0 

0.4251J;®" 



-21.0 

0.716l®;l™ 

-19.0 

^•"^^'^-0.068 



-20.0 

0.678l°;™e 

-18.0 

0.356lJ;™4 


Note. — Columns: (1) (3) (5) UV magnitude. (2) (4) (6) Median ef¬ 
fective radius at the rest-frame optical or UV wavelength. The lower and 
upper limits indicate the 16th and 84th percentiles of the Ve distribution, 
respectively. 

(e.g., Haufiler et al. 2013; Vika et al. 2013, 2014). Vul- 
cani et al. (2014) have reported that n tends to be larger 
in redder bands for blue galaxies due to a bulge compo¬ 
nent with old stellar ages and/or dust attenuation at the 
central region. 

Our results confirm that n values of QGs are signifi¬ 
cantly higher than those of SFGs at z < 2 in the second- 
top right panel of Figure 5. For massive SFGs with 
logM* = 10 — 11 M 0 , n values monotonically increase 
from n^l — 1.5atz~lto'n^2 — 3 at 2 :~ 0 . The 
evolutionary trend of n for the massive SFGs is similar to 
that of the QGs at z ~ 0 — 2, which is consistent with pre¬ 
vious results (see the discussions in Pastrav et al. 2013; 
Naab et al. 2009; van Dokkum et al. 2010). 

A.t 2 < 2 ; < 3, n values of the SFGs at the rest-frame 
optical wavelengths are smaller than those at the rest- 
frame UV wavelengths slightly by An < 0.5, which is 
similar to the results of Vulcani et al. (2014) for local 
objects. 

Interestingly, in Figure 5, we find that typical SFGs 
have a value of n ^ 1 — 1.5 at the wide redshift range 
of 2 : ~ 1 — 6 , albeit with the large scatter of individual 



0 12 3 4 

rY^ [kpc] 


Fig. 6 .— Distribution of for the SFGs and the LBGs at 
z 1 — Q in the bin of L\jy = 0.3 — 1 L*_^. The histograms and 
the curves show the distributions and the best-fit log-normal 
functions, respectively, for the SFGs at 2 = 1 — 2 (green) and 2 — 3 
(light-green) and the LBGs at 2 : 4 (blue), 5 (light-blue), and 6 

(cyan). The y-axis is arbitrary. The histograms and curves are 
slightly shifted along x- and y- axes for clarity. The shifted values 
are Are = —0.25, —0.12, —0.09, —0.04, and 0 kpc for 2 : = 1 — 2 and 
2 : = 2 — 3 SFGs, and 2 : 4, 2 : 5, and 2 : ~ 6 LBGs, respectively. 

Although these choices of the shifts moderately cancel out the trend 
of the re evolution, the re decrease towards high- 2 : is still clearly 
found. 

galaxies. There is a similar claim made by e.g. Morishita 
et al. (2014), but only for 2 : ^ 1 — 3 star-forming galaxies 
(Figure 5). Our results newly suggest that the typical 
Sersic indices of star-forming galaxies are n ^ 1 — 1.5 at 
z~3-6. 

This constant n guarantees that we use a fixed n value 
of 1.5 in the size measurements for LBGs (Section 3). 

5.2. Size-Luminosity Relation 

We investigate the size-luminosity rg-Luv relation and 
its dependence on redshift. Figure 9 and Table 8 repre¬ 
sent the size-luminosity relation at z = 0 — 8 for the SFGs 
and LBGs, where Luv is presented with Muv- We can¬ 
not examine the size-luminosity relation at z ~ 10, be¬ 
cause the number of z ^ 10 LBGs is only three. A large 
area of ^ 910 arcmin^ in the HST fields allows us to de¬ 
rive the re-L\jy relation in a wide range of magnitude, 
—23 < Muv ^ —17 mag even for z ~ 4 LBGs. Figure 
9 shows that re has a negative correlation with Muv at 
0 < 2 < 8 . 

The Te-Luv relation is fitted by 



(6) 
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Fig. 7. — Distribution of Ve in different .Luv bins, 0.12 — 0.3 (left), 0.3 — 1 (middle), and 1 — 10 (right) L\j\/L*_^. Each row displays 
galaxies from z = 0 — 1 (bottom) to 2 : = 8 (top). The red, green, and blue histograms indicate distribution of and for the SFGs, 
and for the LBGs, respectively. The solid curves denote the best-fit log-normal functions for these histograms. The solid and dashed 
arrows present the median and average values of re with the color coding same as the curves. The y-axis is arbitrary. 
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Redshift 


Fig. 8. — Difference of the size evolution results based on the average (crosses), median (filled circles), and modal (open diamonds) values 
of re in the bin of Lyv = 0-3 — The red, cyan, and blue filled symbols indicate r®*’* and for the SFGs and for the LBGs, 

respectively. The error bars for our re are not plotted for clarity, although these errors are included for estimating fiz- The solid, dashed, 
and dot-dashed curves denote the best-fit size evolution for the average, median, and modal re values, respectively, in the linear space. 
Note that these differences of the statistical results are found in the linear space of re, because the re-distributions follow the log-normal 
functions (see the text). The re values for LBGs in the literature are plotted with gray symbols (open circles; Gurtis-Lake et al. 2014, open 
pentagon; Holwerda et al. 2014, open diamonds; Huang et al. 2013, filled triangles; Ono et al. 2013, open triangle; Grazian et al. 2012, 
open inverse-triangles; Oesch et al. 2010, and open squares; Bouwens et al. 2004). 


where tq and a are free parameters. The rg value repre¬ 
sents the effective radius at a luminosity of Lq, which is 
similar to the parameter 7 used in e.g., Newman et al. 
(2012). The a value is the slope of the r^-Luy relation. 
We select Lg to the best-fit Schechter parameter M* at 
z ~ 3 that corresponds to Myv = —21.0, following the 
arguments of Huang et al. (2013). 

The left panel of Figure 10 shows the redshift evo¬ 
lution of rg and a. We parametrize the size growth 
rate by fitting rg with a function of Bz{l + z)^’‘. The 
best-fit function is 6.9(1 -I- 2)“i'20±-o.o4 which do 
not significantly change even with and without the 
results. We also carry out fitting with a function of 
where Bh and Ph are free parameters and 
h{z) = H/Hg = A/Om(l + zy + Ha. Here the fit¬ 
ting of the /i(z)-form functions are conducted, because 
these /i(z)-form functions could be a more realistic phys¬ 
ical treatment as claimed by van der Wei et al. (e.g., 
2014). The fitting results yield the best-fit function of 
5.3kpc that is plotted in the left panel of 
Figure 10. Although we do not use the Ve estimate of 
z ^ 10 for the fitting, the z ^ 10 data point is placed on 
the the extrapolation of the best-fit function. 

The evolution of rg is similar to those of the median r^ 
values that are presented in Figure 8 . Here we plot re as 
a function of redshift in Figure 11, which is the same as 
Figure 8 , but for the median re values of three different 


UV luminosity samples. We fit the functions and find 
that the best-fit are —1.22 ± 0.05, —1.10 ± 0.06, and 
-0.84 ± 0.11 in Luv/Ll^g = 0.12 - 0.3,0.3 - 1, and 
1 — 10, respectively (Table 6). The best-fit f3z values are 
comparable to the one of rg. 

In contrast to the re evolution, there is no significant 
evolution of a (eq. 6) at z = 0 — 8 found in the right 
panel of Figure 10. We calculate the weighted-average 
value of a with our data points over z = 0 — 8, and 
obtain a = 0.27 ± 0.01. Figure 10 compares the a es¬ 
timates of z = 0 — 8 obtained in the previous studies. 
The OL measurements of local spiral and/or disk galaxies 
are comparable to a ^ 0.27 (Shen et al. 2003; Courteau 
et al. 2007; de Jong & Lacey 2000). At z = 0 — 3, van 
der Wei et al. (2014) have revealed that the slopes of 
size-stellar mass relation do not evolve. Adopting eq. 
(2) to calculate Luv from the stellar masses, we obtain 
the rg-Luv relation and evolution similar to our results. 
At z > 4, there are several a measurements reported by 
Curtis-Lake et al. (2014); Huang et al. (2013); Jiang et al. 
(2013); Grazian et al. (2012). However, these data points 
of a are largely scattered (the right panel of Figure 10). 
Nevertheless, our a values fall within the scatter of the 
previous measurements. 

Our results of the re (or rg) evolution and the con¬ 
stant OL suggest that the rg-Luv relation of star-forming 
galaxies is unchanged but with a decreasing offset of re 
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Fig. 9.— Effective radius re and UV magnitude Muv relation at 2 ~ 0 — 8. The top, middle, and bottom panels represent r^^* and 
for the SFGs, and r^^^ for the LBGs, respectively. The redshifts for the relations are labeled at the top of the panels. The red lines denote 
the best-fit power-law functions of oc Tuv for tho re-Muv relations. The best-fit power law for the z ~ 0 — 1 SFGs are plotted in the 
all panels (the dashed lines). The open squares in the bottom panels denote values obtained with the stacked images of LBGs for the 
purpose of the cosmological SB dimming effect evaluation (see Section 5.2). The gray points with error bars indicate the median and 
the 16th and 84th percentiles of the distribution. 
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Fig. 10.— Results of power-law fits for the re-Muv relation in Figure 9. The red, cyan, and blue filled circles indicate estimates of ro 
and a based on and for the SFGs, for the LBGs, respectively. (Left) Effective radius ro at Luv = 1 3 corresponding 

to Muv = —21. The thin dashed, dotted and thick dashed lines show the best-fit curves with (1 -|- z)~^, (1 -|- and (1 -|- , 

respectively. The dot-dashed line indicates the fit of rg oc . The best-fit l3z and I3h values are —1.20 di 0.04 and —0.97 di 0.04, 

respectively. (Right) Slope o; of rg oc Luv as a function of redshift. The thick dashed and thin gray lines denote the weighted-average 
value with a Icr error, a. = 0.27 di 0.01. The open symbols show cx. for the SFGs or the LBGs in the literature (open squares assuming 
Equation 2; van der Wei et al. 2014, open circles; Curtis-Lake et al. 2014, open diamonds; Huang et al. 2013, cross; Jiang et al. 2013, and 
open triangle; Grazian et al. 2012). The gray filled symbols represent the results for local spiral and/or disk galaxies (filled triangle and 
inverse-triangle; n < 2.5 galaxies with r-band magnitudes of Mr < —20.91 and Mr > —20.91, respectively, in Shen et al. 2003, filled circle; 
Courteau et al. 2007, and filled square; de Jong & Lacey 2000). 
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Fig. 11.— Redshift evolution of median rg in different Luv bins, 0.12 — 0.3 (left), 0.3 — 1 (center), and 1 — 10 (right) Luv/T*_ 3 . The 
definitions of the symbols and lines are the same as those in Figure 10. The data points are slightly shifted along x-axis for clarity. 


from 2 ; = 0 to 8. Because the morphological evolution 
trend of star-forming galaxies is simple, our results ben¬ 
efit to studies using Monte-Carlo simulations for lumi¬ 
nosity function determinations that require an assump¬ 
tion of high-z galaxy sizes (e.g., Ishigaki et al. 2014; 
Oesch et al. 2014). Moreover, these morphological evo¬ 
lution trends are important constraints on parameters of 
galaxy formation models. 

Note that there is a possible source of systematics 
given by the cosmological SB dimming effect by which 
we would underestimate Vg (Section 3). To estimate the 
effect of the cosmological SB dimming, we measure 


of z ~ 4 — 8 LBGs with stacked images that accomplish 
the detection limit deeper than the individual images by 
a factor of ~ 20 — 30. The values measured in the 
stacked images roughly reproduce the size-luminosity re¬ 
lation of Figure 9, suggesting that there are no signatures 
of systematics in the values measured by our GALFIT 
profile fitting technique. There is another possibility of 
the cosmological SB dimming effect. If there exist a large 
population of diffuse high-z galaxies that are not identi¬ 
fied in our HST images, we would underestimate the 
values. However, it is unlikely that such a diffuse high-z 
population exists. This is because the luminosity func- 
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Fig. 12.— SFR SD SsFR as a function of redshift. The definitions 
of the symbols are the same as those in Figure 10. The large, 
middle, and small circles denote Ssfr in the Luv bins of 1 — 10, 
0.3 — 1, and 0.12 — 0.3 -Luv/T*^ 3 , respectively. The filled circles 
are the same as in Figure 10. The SFR for the LBGs is corrected 
for dust extinction with two relations of (Bouwens et al. 

2014a) and IRX-/3 (Meurer et al. 1999). The dashed gray lines 
represent the Ssfg evolution calculated with an SFR of 10 Mq/jy 
and the best-fit rp curve in Figure 10. The open symbols are taken 
from the literature on LBGs (triangles; Ono et al. 2013, inverse- 
triangles; Oesch et al. 2010). The error bars denote the 16th and 
84th percentiles of distribution. The data points are slightly shifted 
along x-axis for clarity. 

tions of z ^ 4—6 LBGs derived with HST data agree with 
those obtained by ground-based observations (Beckwith 
et al. 2006) whose PSF’s FWHM is ^ 1" corresponding 
to 3 — 4 kpc in radius at z ~ 4 — 6. In other words, 
at these redshifts, there is no diffuse population with a 
radius up to ^ 3 — 4 kpc that is significantly larger than 
our size measurements of rg < 1 kpc (see, e.g.. Figure 6). 
We therefore conclude that our results of size measure¬ 
ments are not significantly changed by the cosmological 
SB dimming effect. 

5.3. SFR Surface Density 

We examine the redshift evolution of SFR SD Ssfr- 
Figure 12 shows Ssfr as a function of redshift. Figure 12 
is the same as Figure 5, but for all of our galaxies up to 
z = 8 with the binning of Luv values. Figure 12 shows 
that SsFR gradually increases by redshift from z ^ 0 
to 8. This evolutional trend and the Ssfr values are 
consistent with those of z ~ 4 — 8 previously reported by 
(e.g., Oesch et al. 2010; Ono et al. 2013). Our results of 
the SsFR evolution suggests that Ssfr of typical high-z 
galaxies continuously increases from z ~ 0 to 8. 

In Figure 12, we also find that the increase rate per red¬ 
shift becomes small at z > 4 in the regime of log Ssfr ~ 
0.5 — 1 Mq yr“^ kpc“^. We obtain the Ssfr evolution 
curve using the eq. (3) with the inputs of the best-fit 
function rg = 6.9 (i q- 2 :)“i' 20±-0 04 (Section 5.2) and the 
SFR estimated from the Luv value via equation (4). Fig¬ 
ure 12 presents the Ssfr evolution curve. As expected, 
the Esfr evolution curve follows the Ssfr data points. 
In other words, the slow Ssfr evolution at z > 4 is ex¬ 


plained by the simple power-law galaxy size evolution of 
rg = 6.9 (1 -I- 

In Figure 13, we examine the dependence of Esfr on 
SFR and M,. The left and right panels of Figure 13 show 
SsFR as functions of SFR and M,, respectively. For com¬ 
parison, we also plot SDSS galaxies with an exponential 
SB prohle in Lackner & Gunn (2012) and the Milky- 
Way (Kennicutt & Evans 2012). These local galaxies are 
placed in the regime of low Ssfr values. Obviously, Fig¬ 
ure 13 reproduces the result of Figure 12 that Esfr is 
typically higher for high-z galaxies than low-z galaxies. 
In the Esfr-SFR diagram of Figure 13, Esfr positively 
correlates with SFR. This is because the Esfr and SFR 
values are related by eq. (3). The slopes of Esfr-SFR 
relation appear similar at z ^ 2 — 8. On the other hand, 
we find that the Esfr-AT* diagram of Figure 13 shows no 
strong dependence of Esfr on M* (see also, e.g. Wuyts 
et al. 2011). These two diagrams suggest that Esfr in¬ 
creases towards high-z, keeping the similar Esfr-SFR 
and Esfr-AI* relations over z ~ 2 — 8. 

5.4. Size-UV Slope P Relation 

We derive the rg-UV slope /3 relation to investigate the 
dependence of galaxy sizes on stellar population. The /3 
parameter is defined by f\ oc where fx is a galaxy 
spectrum at ~ 1500 — 3000A, which is a coarse indicator 
of the stellar population and extinction of galaxies. A 
small P means a blue spectral shape, suggesting young 
stellar ages, low metallicity, and/or dust extinction. 

For the SFGs, we calculate P via 

g _ 4741700 — ni2800 2 (J) 

^ “ ~2.51ogl700/2800 ~ ^ ^ 

where mi 7 oo and m 28 oo are the total magnitudes at wave¬ 
lengths of 1700 and 2800 A in the rest-frame, respec¬ 
tively. These magnitudes are taken from the catalogue 
of Skelton et al. (2014). For the z ~ 4, 5, and 6 LBGs, 
we derive /3, fitting the function of fx oc to the mag¬ 
nitude sets of *775/8142:850^105'^125 j 2 : 8501^105 Ji25L7i 60, 
and yio 5 Ji25L7i 60, respectively, in the same manner as 
Bouwens et al. (2014a). For the z ^ 7 and 8 LBGs, we 
estimate P using 

P= -2.0 + 4.59(Ji25 - Hieo) (forz - 7), (8) 

P = -2.0 -f 8.68( JL1 i4o - i7i6o) (for z - 8). (9) 

Figure 14 represents the Ve-P relation in the bin of 
Luv/A *^3 = 0.3 — 1. We hnd that Luy-beta relation 
is poorly determined for the z ^ 8 LBGs, due to the 
small statistics, and the z ~ 8 result is not presented. In 
Figure 14, we identify clear trends that smaller galaxies 
have a bluer UV spectral shape at 0 < z < 7. This is 
consistent with the results of z ~ 6 — 8 LBGs reported by 
Kawamata et al. (2014). This Vg-P correlation indicates 
that young and forming galaxies have typically a small 
size. We find a negative correlation between Cg and P 
for the z = 5 — 6 SFGs. The negative-correlation trend 
appears simply due to the small sample, which is not 
statistically significant. 

6. DISGUSSION 

6.1. The Tg Distribution and SHSR: 
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Fig. 13.— SFR SD Ssfr as functions of SFR (left) and stellar mass (right). The small magenta and cyan circles indicate median SgpR 
values at a given SFR or M* tor the SFGs at 2 ~ 2 and 2 ~ 4, respectively, based on ■ The large circles represent the LBGs at 2 ~ 4 
(cyan), 2 ~ 6 (green), and 2 ~ 8 (dark-blue). The dark-blue points denote individual LBGs at 2 = 7 — 8. The SFR tor the LBGs is 
corrected for dust extinction with two relations of Mu-y-/3 (Bouwens et al. 2014a) and IRX-/3 (Meurer et al. 1999). The square represents 
the Milky Way (Kennicutt &; Evans 2012). The gray dots indicate SDSS galaxies with an exponential SB distribution from a catalog of 
Lackner &; Gunn (2012) whose Ssfr calculated from the SFR and values based on u-band magnitudes and single Sersic profile fits, 
respectively. The stellar mass of the SDSS galaxies is taken from Kauffmann et al. (2003); Brinchmann et al. (2004); Salim et al. (2007). 
The dashed lines correspond to constant effective radii of = 0.1,1,10 kpc, from top to bottom. The horizontal lines are the weighted 
average values of Ssfr in each redshift bin. The error bars denote the 16th and 84th percentiles of the distribution. 


Implications for Host DM Halos and Disks 

Here we investigate the properties of the re distribu¬ 
tions in Section 6.1.1, and estimate SHSRs in Section 
6.1.2. Combining these results and theoretical models, 
we discuss the host DM halos and the stellar dynamics 
in Section 6.1.3. 

6.1.1. Log-Normal Distribution of r^ 

In Section 4.2, we find that the r^ distributions of our 
galaxies are well fitted by the log-normal functions in the 
wide-range of redshift, z ~ 0 — 6, and luminosity. 

Figure 15 shows the best-fit ainr^ values as a func¬ 
tion of redshift. Size measurement uncertainties crinre.err 
would broaden the width of the re-distribution. We es¬ 
timate typical (Tinre.err in cach 2 and Luv bin. We cor¬ 
rect CTin re for the size measurement uncertainties through 

(Tlnre “ ((Tlnre.obs (Tlnre.err where (Tlnre,obs i® the 

observed width of the rg distribution. We hnd that CTin re 
values fall in the range of ^ 0.45 — 0.75 with no clear 
evolutional trend at 2 ~ 0 — 6. Our a\nr^ values are 
slightly larger than the estimates for local disks in Shen 
et al. (2003); de Jong & Lacey (2000); Courteau et al. 
(2007) and for late-type galaxies at 2 ^ 0 — 3 in van 
der Wei et al. (2014). These differences would be ex¬ 
plained by the choices of the wavelengths for the galaxy 
size measurements, because these previous studies mea¬ 
sure galaxy sizes in the rest-frame optical wavelength. 
In fact, if we change from the rest-frame UV-luminosity 
to optical wavelength sizes for the size distribution, we 
obtain moderately small din re values. However, there 
still remain the differences of ^ 20 — 30% beyond the 


error bars in Figure 15. These ^ 20 — 30% differences 
are probably explained by the sample and measurement 
technique differences. We also compare the din re esti¬ 
mates of 2 ~ 4 — 5 LBGs given by Huang et al. (2013), 
and find a moderately large difference by a factor of 1.5. 
However, the scatters of our measurements and the sta¬ 
tistical uncertainties of Huang et al. 2013’s estimates are 
too large to conclude the differences. 

6.1.2. SHSR 

We estimate the SHSRs that are defined with the ratio 
of rgjr^yc^ where Tvir the virial radius of a host DM halo. 

The Tvir value is calculated by 

2GMvir 

Avirflm(-2)-ff(z)^ 

where Avir = IStt^ -|- 82a; — 39a;^ and x = flmiz) — 1 
(Bryan & Norman 1998). We obtain the virial mass of 
a DM halo, Mvn, from stellar mass, M*, of individual 
galaxies by using the relation determined by the abun¬ 
dance matching analyses (Behroozi et al. 2010, 2013). 
Figure 16 shows rgfr^i^ as a function of redshift and its 
dependence on Ljjy at 2 ~ 0 — 8. The 2^10 data point 
is omitted due to small statistics. In Figure 16, we find 
that Tgjryyg is ^ 2% for the star-forming galaxies and 
^ 0.5% for the QGs. Interestingly, rgjryn- of the star¬ 
forming galaxies is almost constant with redshift, albeit 
with the large uncertainties at 2 > 5. The no signifi¬ 
cant evolution of rgjry-a is reported by Kawamata et al. 
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Fig. 14.— Relation between effective radius re and UV slope /3 for the SFGs and the LBGs with Lxjv = 0.3 — 1 L*_ 3 . The symbols are 
the same as in Figure 9. The red lines denote the best-fit power-law functions of Ve oc 0^, where c is a free parameter. The 2 ~ 8 relation 
is not shown here due to the poor statistics. 


(2014) based on a compilation of data from the literature 
for star-forming galaxies at z > 2. Our systematic struc¬ 
tural analyses confirm the report of no large evolution 
seamlessly from z ~ 0 with the homogenous data sets 
and the same analysis technique over the wide redshift 
range. Figure 16 also indicates that there is no strong 
dependence of Te/cvir in the wide luminosity range of 
Luv ~ 0.12- 10 L*^ 3 . 

We compare our rg/rvir estimates with those of pre¬ 
vious studies. Because the previous studies choose dif¬ 
ferent statistics for Tg/rvir estimates, we present aver¬ 
age, median, and modal r-g/r-vir for our galaxies with 
0.3 — lFuv/i ^=3 in Figure 17, together with the pre¬ 
vious results. 

For local galaxies, Kravtsov (2013) obtain rg/fyir = 
1.50 ± 0.07% by the fitting of size-luminosity relations. 
This result of z = 0 is consistent with our results at 
a similar redshift of z ^ 0.5 within the Icr uncertainty 
(Figure 17). For high-z galaxies, Kawamata et al. (2014) 
calculate rg/cvir values with the average statistics. In 
Figure 17, the gray symbols of Kawamata et al. 2014’s 
estimates agree with blue crosses of our results. We find 
that the results of ours and the previous studies fall in 
the Tg/cvir range of r-g/cyir = 1.0 — 3.5%, regardless of 
statistics choices. 

Motivated by the no large evolution of rejr^h-, we cal¬ 
culate (rg/ryir) that is a value weighted-averaged 

over z ^ 0 — 8. We obtain (rg/ryir) = 2.76 ± 0.47%, 
1.92 ± 0.09%, and 1.13 ± 0.06% for our Cg/cyir estimates 
of average, median, and modal statistics, respectively. 

The (rg/cyir) value from our average statistics results is 


in good agreement with that of Kawamata et al. (2014), 
3.3 ±0.1%. 

6.1.3. Dark Matter Halo and Stellar Disk 

Summarizing our observational findings for star¬ 
forming galaxies in Sections 6.1.1 and 6.1.2, we identify, 
over cosmic time of z ^ 0 — 6, that the Cg distribution 
is well represented by log-normal distributions, and that 
the standard deviation is ciinrs 0-45 — 0.75, and that 
the SHSR is almost constant, ^ 2%. It is interesting 
to compare these observational results with the theoret¬ 
ical predictions of the spin parameter A distribution of 
host dark halos. DM N-body simulations suggest that A 
follows a log-normal distribution with the standard de¬ 
viation of (Tin A ~ 0.5 — 0.6 (e.g., Barnes & Efstathiou 
1987; Warren et al. 1992; Bullock et al. 2001). The shape 
and the standard deviation of the A distributions are very 
similar to those of Cg. These similarities support an idea 
that galaxy sizes of stellar components would be related 
with the host DM halo kinematics. Our study has ob¬ 
tained this hint of Cg-A relation at the wide range of red¬ 
shift, z ^ 0 — 6, that complements the previous similar 
claim made for z < 3 galaxies (van der Wei et al. 2014). 

If Tg values are really determined by A as indicated by 
the Tg distribution properties, stellar components of the 
high-z star-forming galaxies have dominant rotational 
motions that form stellar disks. In fact, according to 
disk formation models (e.g.. Fall 1983, 2002; Fall & Ef¬ 
stathiou 1980; Mo et al. 1998), gas receives the specihe 
angular momentum from host DM halos through tidal 
interactions which make a constant SHSR similar to the 
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Fig. 15.— Standard deviation of the log-normal r^, distri¬ 

bution (Equation 5) as a function of redshift. The colored symbols 
are the same as in Figure 10, but the large, medium-large, and 
small circles denote fringe 'u the Luv bins of 1 — 10, 0.3 — 1, and 
0.12 — 0.3 LuyIL*_^, respectively. The shaded region indicates 
the width of the A distribution, y, predicted by N-body simu¬ 
lations (e.g., Bullock et al. 2001). The value at 2 > 7 is not 

plotted due to the poor statistics. The gray open symbols show 
frinre for SFGs Or LBGs in the literature (open squares; van der 
Wei et al. 2014, open diamonds; Huang et al. 2013). The gray filled 
symbols represent the results for local spiral or disk galaxies (filled 
inverse-triangle for n < 2.5 galaxies with Mr > —20.91; Shen et al. 
2003, filled circle; Courteau et al. 2007, filled square; de Jong & 
Lacey 2000). 

one found in Section 6.1.2. 

Moreover, in Section 5.1, we find that typical high-z 
star-forming galaxies have a low Sersic index of n ~ 1.5 
at z ~ 0 — 6. The combination of the log-normal Ve 
distribution, the rg-A standard deviation similarity, and 
the low Sersic index suggests a picture that typical high- 
z star-forming galaxies have stellar components similar 
to disks in stellar dynamics and morphology over cosmic 
time of z ^ 0 — 6. 

6.2. Specific Disk Angular Momentum 
Inferred from the Observations and Models 

As we discuss in Section 6, a number of observational 
results suggest that typical high-z star-forming galaxies 
have disk-like stellar components in dynamics and mor¬ 
phology at z ^ 0 — 6. Thus we compare our results with 
the disk formation model of Mo et al. (1998), 

re 1.678 / Jd ^ /r(A, Cvir j Wd, jd) 

T'vii: \/2 y rUd j \//c(Cvir) 

where 1.678 is a coefficient for converting the scale length 
of exponential disk i?d to r^,. The jd (wd) value is a 
angular momentum (mass) ratio of a central disk to a 
host DM halo. The /c(cvir) and /r(A, Cvh-, rud: jd) are 
functions related to halo and baryon concentrations, re¬ 
spectively. The Cvir is the halo concentration factor. The 
full functional forms of /c(cvir) and /r(A, Cvr, rud, jd) are 
found in Mo et al. (1998). The SHSR Ce/rvir with a fixed 


Fig. 16.— Median SHSR, ^e/rvir, as a function of redshift. 
The colored symbols are the same as in Figure 10, but the large, 
medium-large, and small circles denote rg/rvir values in the L^jv 
bins of 1 — 10, 0.3 — 1, and 0.12 — 0.3 -Luv/^s— 3 ? respectively. 
The horizontal dashed line indicates a weighted mean of (re/rvir) 
in the 0.3 — l-Luv /-^*=3 bin. The red squares denote the QGs 
with logM* = 10.5 — 11 Mq. The virial mass of host DM halos is 
derived from the results of Behroozi et al. (2013). 


jd/^d shows little or no dependence on and jd. If we 
use A and Cvir values well constrained by numerical sim¬ 
ulations (e.g., Vitvitska et al. 2002; Davis & Natarajan 
2009; Prada et al. 2012), we can constrain jd/md- 

Figure 17 presents r-g/rvir regions corresponding to 
jd/^T-d = 0.5, 1.0, and 1.5. To determine these re¬ 
gions, we randomly change the A and Cvir values within 
A = 0.038 — 0.045 (Vitvitska et al. 2002; Davis & Natara¬ 
jan 2009) and Cvr ranges at logMvR = 11—13 Mq in 
Figure 12 of Prada et al. (2012), respectively. We also 
assume the conservative range of 0.05 < < 0.1 (e.g.. 

Mo et al. 1998). Substituting these numbers and our re¬ 
sults of Te/rvir (Section 6.1.2) into eq. (11), we obtain 
jd/m^ = 0.7 — 0.8. Note that our estimates of r^jr^ir 
fall in jd/wd ~ 0.5 — 1 at z ~ 0 — 8, regardless of the 
statistical choices (Figure 17). 

This result of jd/nid "^0.5 — 1 indicates that a central 
galaxy acquire more than half of specific angular mo¬ 
mentum from a host DM halo. Our jd/wd values are 
comparable to the estimates with kinematical data for 
nearby disks {jdl'md 0.8; Romanowsky & Fall 2012; 
Fall & Romanowsky 2013). Moreover, Genel et al. (2015) 
predict jd/wd 1 for z ~ 0 late-type galaxies with the 
Illustris simulations (Vogelsberger et al. 2014a,b; Genel 
et al. 2014). These independent studies for z ~ 0 galaxies 
confirm that our estimate of jd/m-d ^ 0.5 — 1 is correct 
at z ~ 0, and suggest that the conclusion of no signifi¬ 
cant evolution of jd/wd over z ^ 0 — 8 would be reliable. 
Genel et al. (2015) have revealed that galactic winds with 
high mass-loading factors (AGN feedback) enhance (sup¬ 
press) jd/nid. This suggests that the no significant evo¬ 
lution of jd/m-d at 0 < z < 8 would place important 
constraints on parameters of galaxy feedback models. 

In Section 6.1.2, we obtain that the SHSR of QGs is 
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Fig. 17.— Comparison between our SHSR, re/rvir? and those of previous studies in a bin of ~ 0.3 — 1 Luv/T^^— 3 - The symbols and lines 
are the same as those in Figure 8, and we include a measurement for local galaxies with the black filled circle (Kravtsov 2013). The re/Tvir 
values of the gray symbols are taken from Kawamata et al. (2014) who compile the results of the literature. The horizontal dashed, solid, 
and dot-dashed lines indicate weighted means of (r’e/?'vir) of average, median, and modal values, respectively. The red, green, and blue 
shaded areas illustrate the regions of Ja/^d = 1-5, 1-0, and 0.5, respectively (see Section 6.2 for details). A typical error bar in our re/^vir 
estimates is shown at 2 : 0.5. 


^ 0 . 5 % that is about four times smaller than the one 
of star-forming galaxies. If we naively assume that QGs 
follow eq. (11) with the one-forth of the specific an¬ 
gular momentum of the star-forming galaxies, we ob¬ 
tain jd/nid ^ 0.1 — 0.25. This value is comparable to 
jd/rrid ~ 0.1 for nearby ellipticals in Fall & Romanowsky 
(2013) and jd/uid ^ 0.3 for 2 ~ 0 early-type galaxies pre¬ 
dicted in Genel et al. (2015). This small specific angular 
momentum of QGs would be explained by the loss of an¬ 
gular momentum via dynamical frictions during merger 
events and/or weak feedback (e.g., Scannapieco et al. 
2008; Zavala et al. 2008). 

6.3. Clumpy Structures of 
High-z Star-Forming Galaxies 

Our study has shown a wide variety of morphologi¬ 
cal measurement results, supplemented by the theoret¬ 
ical models. It should be noted that these results are 
based on the structural analyses for major stellar com¬ 
ponents of the galaxies, because we mask substructures 
such as star-forming clumps (e.g., Guo et al. 2014; Mu- 
rata et al. 2014; Tadaki et al. 2014) in our analyses. The 
signatures of the morphological variety could be emerged 
in dispersions of internal colors and SB profiles in re¬ 
cent structural analyses at z ~ 2 — 3 (e.g., Morishita 
et al. 2015; Boada et al. 2015). Moreover, we find that 
the SFR SD, Ssfr, increases towards high-z in Figures 


12 and 13. This fact suggests that star-forming galax¬ 
ies at high-z would tend to have a high gas mass den¬ 
sity, if we assume the Kennicutt-Schmidt law (Kennicutt 
1998b). The gas-rich disks may enhance formations of 
star-forming clumps through the process of disk instabil¬ 
ities (e.g., Genzel et al. 2011). The detailed analyses and 
results of clumpy stellar sub-components for our galaxy 
samples are presented in the paper II (T. Shibuya in 
preparation). 

7. SUMMARY AND CONCLUSIONS 

We study redshift evolution of rg and the size-relevant 
physical quantities such as Sersic index n, r^ distribu¬ 
tion, rgjr^yc^ and r^ — Tuv relation, using the galaxy 
samples at z = 0 — 10 made with the deep extra- 
galactic legacy data of HST. The HST samples consist 
of 176,152 galaxies with a photo-z at z = 0 — 6 from 
the 3D-HST-I-CANDELS catalogue and 10,454 LBGs at 
z = 4 — 10 selected in CANDELS, HUDF09/12, and par¬ 
allel fields of HFF, which are the largest samples ever 
used for studies of galaxy size evolution in the wide- 
redshift range of z = 0 —10. Our systematic size analyses 
with the large samples allow us to measure galaxy sizes 
by the same technique, and to evaluate the biases of mor¬ 
phological A-correction, statistics choices, and galaxy se¬ 
lection as well as the cosmic SB dimming. Using our 
galaxies at z ~ 2 — 3, we confirm that these biases are 
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small, < 30%, in the statistical sense for star-forming 
galaxies at high- 2 , which do not change our conclusions 
of size evolution. 

Our findings in this study are as follows. 

1. The best-fit Sersic index shows a low value of n ~ 
1.5 for the star-forming galaxies at 2 ^ 0 — 6. The 
low n values indicate that a typical star-forming 
galaxy has a disk-like SB profile. 

2. We derive the Ce-Luv relation for star-forming 
galaxies over the wide-redshift range of 2 = 0 — 8. 
The power-law fitting of rg = ro(Tuv/^o)“ reveals 
that ro values significantly decrease towards high- 2 . 
Similar to the evolution of rg, the average, median, 
and modal re values in the linear space clearly de¬ 
crease from 2 ~ 0 to 2 ~ 6. The Ve values in 
any statistics evolve with Ve (x {1 + 2 )“^ °" 

The slope a of the relation has a constant value 
of a = —0.27 ± 0.01 at 0 < 2 < 8, providing an 
important constraint for galaxy evolution models. 

3. The SFR surface density, Esfr, increases from 
2 ~ 0 to 2 ~ 8, while we find no stellar-mass depen¬ 
dence of Esfr in this redshift range. The increase 
of Esfr suggests that high -2 star-forming galaxies 
would have a gas mass density higher than low- 
2 star-forming galaxies on average, if one assumes 
that the Kennicutt-Schmidt law does not evolve 
significantly by redshift. 

4. We identify a clear positive correlation between Ve 
and (3 for star-forming galaxies at 2 ~ 0 — 7 in the 
luminosity range of 0.3—1.0 L*^^. This is explained 
by a simple picture that galaxies with young stel¬ 
lar ages and/or low metal-|-dust contents typically 
have a small size. 

5. The Te distribution of UV-bright star-forming 
galaxies is well represented by log-normal func¬ 
tions. The standard deviation of the log-normal 
Te distribution, ainr^ is ~ 0.45 — 0.75, and ainr^ 
does not significantly change at 2 ^ 0 — 6. Note 
that the structure formation models predict that 
the distribution of a DM spin parameter A follows 
a log-normal distribution with the A-distribution’s 
standard deviation of ain\ ~ 0.5 — 0.6. The distri¬ 
bution shapes and standard deviations of re and A 
are similar, supporting an idea that galaxy sizes of 
stellar components would be related with the host 
DM halo kinematics. 

6. Combining our stellar re measurements with host 
DM halo radii, Tvir, estimated from the abundance 
matching study of Behroozi et ah, we obtain a 
nearly constant value of rg/rvir = 1.0 — 3.5% at 
0 ^ - 2 : < 6 in any statistical choices of average, 
median, and mode. 


7. The combination of the log-normal rg distribution 
with (Tinre 0.45 — 0.75, and the low Sersic index 
suggests a picture that typical high -2 star-forming 
galaxies have stellar components similar to disks in 
stellar dynamics and morphology over cosmic time 
of 2 ~ 0 — 6. If we assume the disk formation 
model of Mo et al. (1998), our rg/rvir estimates 
indicate that a central galaxy acquires more than a 
half of specific angular momentum from their host 
DM halo, jd/rnd — 0.5 — 1. 


These results are based on galaxies’ major stellar com¬ 
ponents, because we mask galaxy sub-structures such as 
star-forming clumps in our analyses. The detailed anal¬ 
yses and results for the clumpy stellar sub-components 
are presented in the paper II (T. Shibuya in prepara¬ 
tion). We expect that future facilities such as the James- 
Webb Space Telescope, Wide-Field Infrared Survey Tele¬ 
scope, Wide-field Imaging Surveyor for High-redshift tele¬ 
scope, and 30-meter telescopes will obtain deep NIR 
images with a high spatial resolution, PSF FHWM of 
< O'.'l — 072, for a large number of galaxies at 2 ~ 10 and 
beyond. Surveys with these facilities will reveal when 
the size-luminosity relation emerges and whether first 
galaxies fall in the extrapolation of the rg evolution and 
the nearly constant relation of rgjrvir = 1.0 — 3.5% at 
2^0 — 8 and 0 — 6, respectively, that we find in this 
study. 
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